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ABSTRACT 

A variety of current models for gamma-ray bursts (GRBs) suggest a common 
engine - a black hole of several solar masses accreting matter from a disk at a 
rate 0.01 to 10 s~^. Using a numerical model for relativistic disk accretion, 
we have studied steady-state accretion at these high rates. Outside about 
10^ cm, the disk is advection dominated; energy released by dissipation is carried 
in by the optically thick gas and the disk does not cool. Interior to this radius, 
for accretion rates greater than about 0.01 Mq s~^, a global state of balanced 
power comes to exist between neutrino losses, chiefly pair capture on nucleons, 
and dissipation. As a result of these losses, the temperature is reduced, the 
density raised, and the disk scale height reduced compared to the advective 
solution. The sudden onset of neutrino losses (owing to the high temperature 
dependence) and photodisintegration leads to an abrupt thinning of the disk 
that may provide favorable geometry for jet production. The inner disk remains 
optically thin to neutrinos for accretion rates up to about 1 Mq s~^. Energy 
emitted in neutrinos is less, and in the case of low accretion rates, very much 
less, than the maximum efficiency factor for black hole accretion (0.057 for no 
rotation; 0.42 for extreme Kerr rotation) times the accretion rate, Mc^. Neutrino 
temperatures at the last stable orbit range from 2 MeV (no rotation, slow 
accretion) to 13 MeV (Kerr geometry, rapid accretion) and the density from 10^ 
to 10^^ g cm~^. The efficiency for producing a pair fireball along the rotational 
axis by neutrino annihilation is calculated and found to be highly variable and 
very sensitive to the accretion rate. For some of the higher accretion rates 
studied, it can be several per cent or more; for accretion rates less than 0.05 
s~^, it is essentially zero. The efficiency of the Blandford-Znajek mechanism in 
extracting rotational energy from the black hole is also estimated. In light of 
these results, the viability of various gamma-ray burst models is discussed and 



^Max Planck Institut fiir Astrophysik, Garching, Germany, 85740; popham@mpa-garching.mpg.de 

^UCO/Lick Observatory, University of Cahfornia, Santa Gruz, GA 95064; 
woosley@ucoHck.org, cfryer@ucolick.org 



- 2 - 



the sensitivity of the results to disk viscosity, black hole rotation rate, and black 
hole mass explored. A diverse range of GRB energies seems unavoidable and 
neutrino annihilation in hyper-accreting black hole systems can explain bursts 
up to 10^^ erg. Larger energies may be inferred for beaming systems. 

Subject headings: black holes - accretion disks - gamma-ray bursts 

1. INTRODUCTION 

Historically the study of black hole and neutron star accretion has been motivated by 
the study of active galactic nuclei and accreting (stellar mass) x-ray sources. In both cases, 
accretion is Eddington-hmited to small rates Ms 10~^(M/ Mq) yr~^. For an interesting 
range of moderate accretion rates, the energy released by viscous dissipation can be radiated 
away efficiently and consequently the disk is thin. 

For very low or very high accretion rates, however, the cooling efficiency of the disk 
changes and can drop to the point where most of the dissipated energy is not radiated, but 
carried into the hole with the accreting gas in an "advection-dominated accretion flow" 
(ADAF). For low accretion rates, this occurs because the radiative processes become very 
inefficient at low density. The dissipated energy stays in the gas, which becomes very 
hot, and this makes the disk geometrically thick. Recently, ADAFs corresponding to low 
accretion rates and optically thin disks have been studied in some detail (see Narayan, 
Mahadevan, & Quataert 1998 for a review). They have been used to model both accreting 
stellar-mass black holes and the supermassive black hole sources believed to lie at galactic 
centers. 

Advection-dominated flows also occur for high accretion rates if the disk is very 
optically thick. Radiation does not diffuse out in a viscous time scale and so again is carried 
into the hole (Katz 1977; Begelman 1978). The solutions resemble, mathematically, the 
ADAFs at low density since in both cases the gas cannot cool. 

Eventually, however, at extremely high accretion rates, the nature of ADAFs can 
change again due to the onset of efficient cooling by neutrinos. For neutron stars, the 
accretion rate where neutrino losses begin to dominate photons in the energy loss budget 
is about 10-'' M© yr-^ (Chevalier 1993; Fryer, Benz, & Herant 1996). For black holes, the 
limit is much higher because energy can disappear into the hole. As we shall see, it is only 
for accretion rates in excess of about 0.01 per second, that neutrino radiation can cool 
the gas accreting into a black hole on a viscous time scale. We refer to such accretion as 
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"neutrino dominated" . The ADAF solutions one obtains neglecting neutrino emission are 
very different. 

Besides wanting to explore new accretion physics, we are motivated by a desire to 
understand GRBs. Most leading models for GRBs (§2) have in common an engine based 
upon a "hyper-accreting black hole" , a black hole of ~ 2 — 10 accreting mass at a rate 
sufficient to consume from 0.1 to several solar masses within, at most, a GRB time scale 
(on average 20 s for the long complex class of bursts (Fishman & Meegan 1995)). 

For the steady state disks we study, the angular momentum in the initial system 
defines the radius where the disk forms and the viscous time scale at that radius gives the 
approximate duration of the event (though longer events occur if the disk is continually fed 
by an external source such as the collapse of a star). One conclusion of our paper will be 
a range of time scales and accretion rates (viscosity dependent) that should characterize 
various GRB models. We also calculate the variation of temperature, density, disk scale 
height, radial drift velocity and luminosity with radius in the disk. From these we are 
able to infer, albeit approximately, the efficiency for jet production based on neutrino 
annihilation (see also Jaroszyhski 1996) as well as the energy density in the disk and the 
number of orbits experienced by an accreting blob. The latter has implications for whether 
the large magnetic fields that have been invoked for the electromagnetic acceleration of 
jets can accumulate (Blandford & Znajek 1977; Narayan et al. 1992; Hartmann & Woosley 
1995; Meszaros & Rees 1997; Paczyhski 1998). 

Our numerical model is discussed in §3. We solve the shm disk equations in the 
Kerr metric for steady accretion. This differs from some previous work which treated the 
disk as a hydrostatic torus, and approximated the effects of accretion by a sequence of 
tori of decreasing mass (Jaroszynski 1996; Daigne & Mochkovitch 1997). The solutions 
are characterized by four physical parameters: the black hole mass, accretion rate, 
black hole spin parameter a, and disk viscosity parameter a. We include, in addition 
to the emission from pair annihilation, neutrino losses from the electron-positron pair 
capture on nucleons (which we find dominates the cooling), and photodisintegration. The 
equation of state, while simple, contains the effects of radiation, relativistic pairs, and 
degeneracy. Our semi-analytic solutions are thus sufficiently realistic to make qualitative 
predictions regarding the viability of various models for GRBs and to guide the complex 
multi-dimensional hydrodynamical studies that will follow (e.g., Eberl, Ruffert, & Janka 
1998; MacFadyen & Woosley 1998; Fryer & Woosley 1998; Fryer et al. 1998). 

We discuss the results of our calculations in §4 and describe a typical solution. We 
examine solutions for a variety of accretion rates, disk viscosity parameters, and black 
hole spin rates. For each case, we derive neutrino luminosities, which turn out to depend 
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strongly on all of these parameters. Approximations to the accretion rate are given in §5 
for various evolutionary scenarios. 

In Section 6 we calculate the efficiencies of gamma-ray production by neutrino 
annihilation for the various models. We also estimate the efficiencies for energy extraction 
by the Blandford-Znajek process. The final section summarizes our conclusions and 
discusses the imphcations of our hyper-accreting black hole disk solutions for GRB models. 

2. EVOLUTIONARY PATHS 

There are many ways in which a black hole may come to experience rapid transient 
accretion. We expect that each occurs to an appreciable extent in nature; indeed given 
current uncertainties in such quantities as the neutron star kick velocity, common envelope 
evolution and the sizes of stripped helium cores during the late stages of evolution, it is 
uncertain just which predominates (Fryer, Woosley, & Hartmann 1998). It is likely, though, 
that GRBs are a diverse set of phenomena, all having at their heart a hyper-accreting black 
hole, but differing in accretion rate, accretion mass, rotation rate of the black hole, and 
surroundings (Table 1). 

2.1. Merging Neutron Steirs 

The oldest of the binary compact object models for cosmological GRBs, merging 
neutron stars, were first mentioned as possible sources by Paczynski (1986) and Eichler et 
al. (1989). Progenitor double neutron star systems have long been known to exist (Hulsc 
& Taylor 1975) and several evolutionary paths for their formation have been proposed (e.g. 
Srinivasan & van den Heuvel 1982; Brown 1995). For a time it was believed that the burst 
might originate during the merger, but this proved too inefficient (Janka & Ruffert, 1996). 
However, depending upon the equation of state, a result of the merger may be a black hole 
of about 2.5 solar masses surrounded by a disk of approximately 0.1 — 0.2 M0 (Ruffert & 
Janka 1998). 

That accretion from this disk might give rise to GRBs was suggested by Narayan et 
al. (1992). The most detailed calculations to date are by Davies et al. (1994), Ruffert, 
Janka, & Schafer (1995), Janka & Ruffert (1996), and Ruffert & Janka (1998). In particular 
Ruffert & Janka (1998) estimate the mass of the disk to be 0.1 - 0.2 M^, the accretion 
rate about 1 s^^, the neutrino luminosity about 10^^ erg s^^, the energy deposition 
along the polar axes by neutrino annihilation about 5 x 10^° erg s~^, and thus the total 
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energy available to make the burst about 5 x 10^^ erg. Simple equations derived in §5 and 
efficiencies calculated in §6 agree well with these numbers. 

2.2. Neutron Stars Merging With Black Holes 

Paczynski (1991) and Narayan et al. (1992) also suggested that the merger of a neutron 
star with a pre-existing black hole of several solar masses might produce GRBs. In the 
dominant formation scenario of binaries consisting of a black hole and neutron star, the 
black hole is formed via hypercritical accretion during a common envelope phase (Bethe 
& Brown 1998). The resultant black hole has a very low mass (~ 3Mq), very similar to 
the double neutron star systems. Here the situation is somewhat improved because of the 
larger mass of the disk that forms, about 0.5 M© (Eberl 1998; Eberl, Ruffert, & Janka 1998; 
Kluzniak & Lee 1998) and the fact that one already has a black hole to start with. In most 
ways, though, this is just a more energetic analogue of the merging neutron star aftermath. 
Accretion rates and neutrino luminosities are higher because of the larger disk mass. 
The neutrino annihilation efficiency is consequently increased to a few per cent. Eberl et 
al. (1998) find disk masses of 0.5 M^, accretion rates highly variable between 1 and 10 
s~^, neutrino luminosities from 10^^ to 10^^ erg s^^, and efficiencies for neutrino annihilation 
of a few per cent. The total energy available for the burst is about 10^^ erg. Depending 
upon beaming, this may still be inadequate to explain GRB 971214 (Ramaprakash et 
al. 1998) without invoking energy sources other than neutrino annihilation, but it may 
account well for the short hard bursts in the BATSE sample (Fishman & Meegan 1995). A 
concern both for this model and the merging neutron star model is the distance the system 
may travel before merging, perhaps too far to be associated with star formation. However, 
because the trigger threshold on Beppo-Sax exceeds a few seconds, we have no information 
on the optical counterparts to short GRBs. 

2.3. CoUapsatrs 

The mechanism whereby the collapsing iron core of a massive star gives rise to an 
outgoing shock and makes a supernova remains controversial (e.g., Mezzacappa et al. 1997; 
Fryer 1998). Surely many massive stars, in their deaths, produce supernovae by forming 
neutron stars in spherically symmetric explosions, but perhaps not all do. Especially as the 

mass of the star becomes more than ~ 30 Mq on the main sequence (pre-supernova helium 
core masses bigger than 12 M^), the accretion rate of the mantle onto the proto-neutron 
star becomes so great, ~ 1 Mq s~^, and the mass of the iron core so large 2.0 Mq), that 
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a black hole may form before a neutrino powered shock is successfully launched. Woosley 
(1993, 1996) suggested that the outcome of such collapses would be gamma-ray bursts. 
This model is currently being explored in two-dimensional calculations by MacFadycn & 
Woosley (1998) and a very similar model, named the "hypernova", has been discussed by 
Paczyhski (1997). An explosion of this sort may have been associated with Supernova 
1998bw (Woosley, Eastman, & Schmidt 1998). 

Early on, material accretes rapidly though the poles and the angular momentum is 
so low that accretion through the equator is also almost unimpeded. The hole rapidly 
grows to over 3 M^. But as the angular momentum of the accreting material rises to 0.5 
- 1 xlO^^ cm^ s^^, a portion of the remainder accretes through a disk. Two-dimensional 
calculations by MacFadyen & Woosley (1998) show that, for a 14 helium star after 
about 10 seconds, the polar regions are evacuated while disk accretion continues at over 
0.1 Mq for about 10 s (with a subsequent long, slow decline). During most of this 
accretion the Kerr parameter of the black hole is over 0.9. 

2.4. White Dwarfs Plus Black Holes 

White dwarfs may also merge with black holes at an appreciable rate (Fryer et al. 
1998). Except for the smaller secondaries, formation scenarios for binaries consisting of a 
white dwarf and black hole parallel those of neutron stars and black holes. Fryer, Woosley, 
& Hartmann (1998) find that the white dwarfs formed in these binaries will preferentially 
have large masses (~ IMq). These massive white dwarfs are tidally disrupted at a few times 
10^ cm and form an accretion disk around the black hole. The accretion into the hole occurs 
by way of a disk at a rate of about 0.01 — 0.07 Mq s^^ and lasts about 10 — 70 s. The lower 
accretion rate and longer timescale corresponds to low-mass white dwarfs (Mwd = O.7M0), 
whereas the high mass white dwarfs which dominate these binaries have higher accretion 
rates, but lower accretion timescales. Again the black hole is spun up by the accretion (if it 
were not already rotating rapidly) to spin parameters of a ~ 0.5. 

2.5. Black Holes and Common Envelope Evolution 

In many of the formation scenarios for compact binaries, a common envelope phase is 
required to create the short periods. The more massive component first makes a black hole 
or neutron star and is later enveloped as the less massive star becomes a giant. A fraction of 
these systems will not eject the hydrogen common envelope; instead, the rapidly accreting 
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compact primary will merge with the secondary's helium core, accreting quickly enough to 
become a black hole if it is not already. The accretion of the helium core onto the black 
hole has been proposed as a GRB progenitor (Fryer & Woosley 1998). The helium core is 
disrupted by the black hole into an accretion disk with a radius equal to a fraction of the 
initial core, ~ 10'' — 10^*^ cm. The accretion rate, especially along the poles, can be quite 
large, initially perhaps 1 Mq s~^, but again the disk accretion rate is expected to be smaller 
(and last longer). 

Detailed calculations do not exist, but the angular momentum is larger than in the 
collapsar model, comparable in fact to the white dwarf- black hole merger. The disk will 
thus form at a large radius, comparable to that of the helium core. A crude estimate is 1 - 
10 xlO^ cm. At this distance the viscous time scale will range from a fraction of a minute 
to 10 minutes (§5; Table 1) and the accretion rate spans a comparably large range. Again 
the black hole may accrete an amount comparable to its mass and be spun up to high Kerr 
parameters. 

3. COMPUTATIONAL PROCEDURE 

3.1. The Physical Model 

The evolutionary scenarios described above all produce rapid disk accretion 
(0.01 — 10 Mq s~^) onto a black hole of a few solar masses. For present purposes, we are 
interested in the gross properties of these disks rather than their detailed time-dcpcndcnt 
behavior; so, for simplicity, we use a steady-state disk model to study this transient event. 
This should be a reasonable approximation, since the viscous time scales in the inner disk, 
which should have the highest temperature and neutrino flux, are much shorter than those 
in the outer disk, on which M should vary. 

As discussed in §1, we expect that the disk will be unable to cool efficiently via photons, 
and so the energy dissipated in the disk will be advected inward. When the temperature 
and density are sufficiently high, neutrino cooling will become efficient. Thus our model for 
the disk includes advection, but also allows for cooling which varies substantially with disk 
radius, unlike most previous models of ADAFs. Also, because the inner regions of the disk 
and the spin of the black hole may be important, we solve for the structure of the disk in 
the Kerr metric. 

Our disk model is based on the advection-dominated disk model of Gammie & Popham 
(1998, hereafter GP) and Popham & Gammie (1998, hereafter PG). This model solves the 
"slim disk" equations in Kerr geometry. The slim disk equations are a more sophisticated 
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version of the standard thin disk equations which include radial pressure gradients and 
radial energy transfer. For a detailed description of these equations and the form they take 
in the Kerr metric, we refer the reader to GP; however, we review them here for the reader's 
convenience. 

The units are such that G — M — c— 1, where M is the mass of the black hole. Later, 
when we present our results, we plot the variables scaled to cgs units. 

The continuity equation is: 

AnrpHV (y^) = (3-1) 

Here r is the Boycr-Lindquist radius, p is the rest mass density, H is the disk thickness, V 
is the radial velocity measured in a co-rotating frame, V = 1 — 2/r + a? /r'^ is a relativistic 
correction factor, and M is the rest mass accretion rate. 

The gas energy equation is 

,vf^V'7^-iit) (3-2) 



Here u is the internal energy, p is the pressure, $ is the dissipation function, and qtot is 
the total coohng rate, as described below. The parameter / measures the importance of 
cooling; if / = 1, all the dissipated energy is advected with the flow, while if / = 0, cooling 
and dissipation occur at the same rate. Most ADAF models have taken / to be constant 
with radius (usually / = 1), but here wc compTitc / directly from the local dissipation and 
cooling rates. Note that if the cooling rate exceeds the dissipation rate, one can have / < 0. 



The radial momentum equation is 



^ fr-^t (3-3) 



1 — V^ dr ^ prj dr^ 
where 

The fr term combines the effects of gravity and rotation, where A = 1 + a'^/r'^ + 2a^/r^ 
and 7^ = 1 /^(l — V'^)/{r'^A), fl = u'^/u* is the angular velocity, and Q± = ±(r^/^ ± a)~^. 
The radial acceleration, on the left-hand side, is given by the difference between fj. and the 
pressure gradient force, where rj is the relativistic enthalpy rj = {p + p + u) / p. 

The angular momentum conservation equation is 

Mir] - AnHrV^ = Mj. (3-5) 
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Here I is the specific angular momentum of the accreting gas, j — const, is the angular 
momentum accretion rate per unit rest mass accreted. The remaining term gives the viscous 
angular momentum transport rate, where is the viscous stress. The expression for is 
rather lengthy, the calculation of it even more so, and the reader is referred to §4 of GP for 
a full discussion. 

The equation of vertical mechanical equilibrium is 

= (3-6) 

where Vz is an effective vertical frequency. We adopt the expression derived by Abramowicz, 
Lanza, & Percival (1997) for u^: 

= ^4 (3-7) 

where S = —Ut is the "energy at infinity" , which is conserved along geodesies. 

One of the radial energy transfer terms included is the advection of entropy. This 
term can become important in situations where the cooling of the disk is very inefficient, 
so that the viscous time scale is shorter than the cooling time scale. This can occur when 
the optical thickness of the disk, as measured from the midplanc to the surface, is either 
very low or very high. In models of advection-dominated disks, a large fraction fadv of the 
energy dissipated in the disk is stored as entropy and advected inward with the accreting 
gas. GP and PG assumed that fadv was constant with radius, and usually took fadv = 1- 

For the extremely high accretion rates that we arc considering, the optical thickness 
of the disk is very large. Thus the effective temperature is much lower than the midplane 
temperature, and the cooling time is long, so it appears likely that the disk is advection- 
dominated. Highly optically thick disks where advection is important have been modeled 
for the case of FU Orionis objects (Popham et al. 1996), but even there the accretion 
rates (~ 10"^ M© yr'^) are many orders of magnitude smaller than those considered 
in this paper. If radiation were the only source of cooling available, the disk would be 
advection-dominated; however, eventually, the disk reaches a temperature where neutrino 
losses become significant. 



3.2. Neutrino and Photodisintegration Losses 



Cooling by neutrino emission is important in regions of the disk where the temperature 
is sufficiently high. Because of the steep temperature and density dependence of the 



-10- 



neutrino emissivity, we expect the onset of neutrino cooling to occur fairly abruptly at 
some transition radius in the disk. Outside this radius, cooling is inefficient and the disk 
is advection-dominated; inside, neutrinos cool the disk efficiently. Because of this, it is no 
longer feasible to use an fadv which is constant with R. 

We include two types of neutrino losses. The first is neutrino emission due to pair 
annihilation. These are computed from the results of Itoh et al. (1989,1990). A rough 
approximation to these results is given by the simple expression 

q^^p ~ 5 X lO^^Ti*! erg cm-^ s"\ (3-8) 

which can be used to compare pair cooling to the other cooling terms. The second type of 
neutrino cooling we include is due to capture of pairs on nuclei. This is computed according 
to 

q,M = 9.0 X lO^VioTi^i^n^e erg cm'^ s'S (3-9) 

where pio = p/lO^'^ g cm^'^, Tn = T/IO^^ K, and Xnuc is the mass fraction of nucleons 
(Qian & Woosley 1996). To avoid confusing this pair capture loss rate with that from pair 
annihilation, we refer to the former as the "URCA" cooling. Note that URCA cooling will 
dominate when 1.8pio > Tfi. 

Xnuc is zero in the outer disk, but photodisintegration breaks down alpha particles into 
neutrons and protons once T reaches about 10^° K. Xnuc is given by 

X^uc = 30.97pro'^'Ty exp(-6.096/rio) (3-10) 

where this expression gives Xnuc < 1, and Xnuc = 1 elsewhere. The photodisintegration 
process cools the gas according to 

dX 

qphot = lO^'^PioV— erg cm-^ s'^ (3-11) 



For accretion rates greater than ~ 1 s ^, the disk begins to be optically thick to 
its own neutrino emission. We limit the actual neutrino emission to the blackbody limit 
(Mayle 1985): 

L.^ = IctbT' (3-12) 

where is the Stephan-Boltzmann constant. This limit corresponds to the Fermi-Dirac 
blackbody luminosity emitted from the neutrinosphere and is the same for each neutrino 
flavor = ^c-i^h-i^t)- This estimate gives rough agreement with the results of Eberl et 
al. (1998) and allows us to extend our results to accretion rates as high as 10 Mq s~^. 
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3.3. Equation of State 

GP and PG assumed a simple gas pressure equation of state. In the optically thin 
disks they were considering, radiation pressure is unimportant, and the densities are low 
enough that degeneracy pressure was negligible. 

In the present problem, gas, radiation, and degeneracy pressure are all significant. 
Accordingly we write the pressure as 

P = pi?T&±^^ + -aT' + k(^Y\ (3-13) 

4 12 \llej 

where R is the gas constant, K = (27r/ic/3)(3/87rm„)^/^ = 1.24 x 10^^, where m„ is the 
nucleon mass, and /^.g is the mass per electron, which we take to be 2. The three terms 
represent gas pressure from nucleons, pressure from radiation and relativistic electrons 
and positrons, and relativistic degeneracy pressure from electrons, respectively. The 
corresponding expression for the internal energy is 

2 4 4 p 

We have assumed that the gas is pure helium before photo disintegration, hence the factor 
(1 -|- 3X„„c)/4 to account for the change in gas pressure. We have compared the results of 
this expression to those produced by a full EOS routine by Blinnikov, Dunina-Barkovskaya, 
& Nadyozhin (1996), and we find that the total pressure agrees to within 10% for 
10^ < p < W\ 10^ < T < 10^° ^ 



3.4. Solution Method 

We set the outer edge of the disk at 10^ Schwarzschild radii, or r = 2 x 10^ GM/c^, and 
the inner edge just outside the event horizon at r = (1 -|- ^/l — a?) GMjc? ^ where T> — 

At the outer edge we impose two boundary conditions: Q and Cg must equal their 
values in the self-similar advection-dominated solution of Narayan & Yi (1994). Two other 
conditions on the flow are provided by the requirement that the flow pass smoothly through 
two critical points. The flrst is the sonic point r^, where \y\ ~ c^. The second is the 
"viscous point" associated with the flnite propagation speed of viscous effects, where 
\y\ ~ c,^. Associated with each critical point are two conditions that must be satisfled for a 
smooth flow, as well as one degree of freedom, the location of the critical point itself. 
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The final boundary condition normafizes the density (the density appears in the basic 
equations only in the form dlnp/dr). We now have all the boundary conditions required to 
solve the four first-order ordinary differential equations for V,l,p, and T, and to find the 
eigenvalue j. 

The equations are solved using a relaxation method, as described by GP and PG. 

4. RESULTS 

Our solutions are characterized by four parameters: the black hole spin a, the viscosity 
parameter a, the mass accretion rate M, and the black hole mass M. The values of these 
are determined largely by the evolutionary path which provides the mass to the disk. 
The black hole mass and spin depend in part on the mass and angular momentum of the 
progenitor star which originally collapsed to produce the black hole, but both M and a can 
be increased substantially by the rapid accretion we are modeling here. The mass accretion 
rate will depend largely on the evolutionary path which leads to the rapid accretion. 
The viscosity parameter a is poorly known, particularly under the kinds of temperatures 
and densities which are reached in these disks. Since the accretion rate is determined by 
the mass of the disk and the viscous time scale on which it accretes, there should be a 
connection between M and a, as discussed in §5. 

4.1. A Typical Solution 

In order to describe the general features of our solutions, we begin by focusing on a 
"standard model" with a = 0, a = 0.1, M = 0.1 M© s-^ and M = 3 Mq. This model 
is shown in Figures 1 and 2, and repeated in some other figures for comparison purposes. 
In Figure 1, we compare this standard model to two other solutions. The first is a pure 
advection-dominated solution calculated using the same code, but with all cooling terms 
turned off, and including only gas pressure. This solution is identical to the one shown 
by GP and PG, scaled to physical units using the appropriate values of M and M. The 
second is an analytic solution for a thin neutrino-cooled disk described in §5, which assumes 
Keplerian rotation, and that pair capture dominates the cooling and gas pressure dominates 
the pressure. Note that this analytic solution assumes a Newtonian gravitational potential. 
These solutions represent the two extremes of no coohng and highly efficient cooling which 
balances the viscous dissipation at each radius. Thus, we expect that the characteristics of 
our neutrino-cooled disk will fall between those of these two solutions. 
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At large radii, the accretion flow is simply an advection-dominated flow. Figure 1 shows 
that the disk is quite similar to the pure advection-dominated flow at logr = 9. Densities 
and temperatures are too small for neutrino cooling to be significant, while optical depths 
are too large for significant photon cooling. The approximate surface density of the disk 
is S = pH ~ 10^^ g cm~^ at r = 10^ cm, so the disk should be extremely optically thick 
to photons, but optically thin to neutrinos. Figure 2 shows that / = 1 at large radii, and 
equivalently that the cooling time scale is much longer than the accretion time scale, so the 
dissipated energy is advected inward before it can be radiated away. The disk is thick, with 
H R, and has substantially sub-Keplerian rotation, reflecting the importance of pressure 
gradients. 

The first significant cooling occurs in the region from logr = 7.5 — 8.2, where the nuclei 
photodisintegrate into nucleons. The nuclei absorb 10^^ erg g~^; this energy loss is larger 
than the energy added by viscous dissipation, and so the entropy of the gas decreases in this 
zone. This is reflected in the negative value of / in this region (Fig. 2). The disk becomes 
thinner, with H ^ 10 cm at r = 10^-^ cm, and the rotation is much closer to Keplerian. 
Fig. 1 shows that the disk becomes more like the thin neutrino-cooled disk in this zone. 
Fig. 2 shows that the radiation pressure, which was the largest pressure term at logr = 8.2, 
drops in importance as the disk becomes denser and more nucleons are produced, and gas 
pressure, which had dominated in the outer disk, is again the largest term at logr = 7.5. 

At logr ^ 7.5, all of the nuclei have photodisintegrated, and the cooling due to 
photodisintegration shuts off rather abruptly, as shown by the sudden increase in / and 
in the cooling time (Fig. 2). Neutrino cooling, predominantly by pair capture, begins to 
become important here, but at logr — 7.5 it is still well below the dissipation rate. As 
the gas continues to fall inward, the increasing temperature and density produce a rapid 
increase in the neutrino cooling rate, which gradually approaches the dissipation rate. This 
can be seen in the slow decrease in / from one (advection-dominated) toward zero inside 
logr = 7.5. Radiation pressure again becomes the dominant term in the inner disk. Note 
that in the innermost portion of the disk, the dissipation rate briefly becomes negative 
where the stress and shear rate change sign due to the use of a causal viscosity prescription 
(see GP for more details); this produces the feature in / at logr = 6 — 6.4. Near the 
horizon, V is approaching c and Q drops toward uj = 2a/ {r^ + o?r + 2o?) = 0. 

The neutrino luminosity of this solution is L^, = 3.35 x 10^^ erg s~^, of which nearly 
90% is from pair capture (see Table 2 for luminosities). This is less than 2% of Mc^, a 
somewhat lower efficiency than the standard thin disk figure of 5.7%, due largely to the 
fact that much of the dissipated energy goes into the entropy of the gas and is advected 
into the black hole. Note that 2 x 10^^ erg s~^ goes into photodisintegration of the nuclei. 
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The neutrino emission from the disk comes predominantly from between log r ~ 6 — 7. The 
disk is thin to neutrinos; S ~ lO-*^^ g cm~^ in the inner disk, while the opacity to neutrinos 
is about 10~^^ cm^ g~^ giving an optical depth to neutrinos of ~ 0.01. 

4.2. Solutions for Other M, a, a, and M 

Figure 3 shows solutions for M = 0.01,0.1, 1, 10 Mq s~^, all with a = 0, a = 0.1, 
and M — 3 M©. The changes in M produce substantial changes in p and in the neutrino 
luminosity, but the disk height H and radial velocity V show only small changes, and the 
temperature increases slowly with increasing M. The density and neutrino luminosity 
both show a strong dependence on M; this is because the increased density increases the 
neutrino luminosity, which makes the disk thinner and denser. The 0.01 Mq s^^ solution 
has a very low neutrino luminosity Lj, ~ 1.5 x 10^^ erg s~^, while the 1 Mq s^^ solution 
has Lj, ~ 8.5 x 10^^ erg s~^, a factor of more than 5000 larger. At low accretion rates, 
the density and temperature are simply too low to permit effective cooling by neutrino 
emission, and nearly all the dissipated energy, apart from that lost to photodisintegration, is 
advected into the hole. At higher accretion rates, the high temperatures and densities allow 
the dissipated energy to be radiated away efficiently. In fact, both the 1 and 10 Mq s"^ 
solutions have neutrino luminosities in excess of the dissipation rate, because the neutrino 
emission radiates away some of the gas entropy in addition to the dissipated energy. This 
can be seen in the variations in / in Fig. 3; the 0.01 Mq s~^ solution has / — 1 inside 
of the photodisintegration zone, while for 0.1 Mq s~^ / drops gradually, reflecting the 
increasing fraction of the dissipated energy which is radiated away by neutrinos. The higher 
M solutions both quickly reach in this zone, demonstrating that neutrino emission is 
radiating away more energy than the disk is dissipating. 

Figures 4 and 5 illustrate the effects of changes in a. Each shows solutions with 
a = 0.1, 0.03 and 0.01; in Fig. 4 M = 0.1 Mq/ s, while Fig. 5 M = 0.01 Mq s-^ All 
have a = 0, M = 3 Mq. In both cases, decreasing a produces a solution with lower radial 
velocity and higher density. This results in more efficient neutrino cooling, with higher 
neutrino luminosities for the same M. In the 0.01 Mq s~^ solutions (Fig. 5), the change in 
a produces a dramatic change in Lj,, since the a = 0.1 solution is very inefficient. The more 
efficient cooling of the a = 0.01 solutions occurs in the inner disk, so the temperature and 
disk height fall below the a = 0.1 solutions there. 

One interesting feature of the a — 0.01 solutions is that they show maxima in the 
density, temperature, and pressure at around logr = 6.6. These are similar to the maxima 
seen by PG at a similar position in the disk for a solution with a = 0.001 and a solution 
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with / = 0.01. Solutions with small values for a and/or / have relatively small radial 
velocities outside of this radius, and build up high densities and pressures there. Inside of 
the pressure maxima, the gas falls in toward the hole, with large radial acceleration. This 
produces the drop in temperature and density, and therefore in neutrino cooling, seen in 
Figs. 4 and 5. This effect can also be seen, to a lesser extent, in the a = 0.03 solutions and 
in the high-M solutions in Fig. 3. 

Figure 6 shows the effects of increasing the black hole spin a for solutions with 
a = 0, 0.5, 0.95. Changes in a have little effect at large radius. Near the hole, the high-a 
solutions have higher densities and stronger neutrino cooling. The horizon also moves 
to smaller radius as a increases. In the a = 0.5 solution, the neutrino luminosity due to 
annihilation rises dramatically from the a — solution. This results from a small increase 
in temperature, since the annihilation luminosity varies approximately as T^. It appears 
that here cooling due to annihilation is acting as a thermostat which prevents T from rising 
any further. In the a — 0.95 solution, most of the increase in neutrino luminosity comes 
from pair capture, due to the dramatic increase in density in the inner part of the disk. 

Figure 7 shows solutions with a = 0.5 for four accretion rates: M = 0.01, 0.1, 1.0, and 
10 M0 s~^. This can be compared to Figure 3, where a — solutions were shown for the 
same accretion rates. The a = 0.5 solutions have higher temperatures and densities, and 
produce higher neutrino luminosities. These solutions might be appropriate to the later 
stages of a burst event, when a substantial amount of mass and angular momentum has 
been added to the black hole. 

Figure 8 shows how a increases as a black hole accretes mass and angular momentum 
from a thin disk; the rates at which the black hole mass and angular momentum increase 
are given by Bardeen (1970). If the black hole is spun up by accretion from an ADAF, it 
will spin up more slowly, and will reach an equilibrium value of a where it neither spins up 
nor spins down; see PG for details. Spinup due to accretion from a neutrino-cooled disk 
will be intermediate between the thin disk and the ADAF cases. 

We have also calculated solutions with larger values of the black hole mass. These 
include a 6 M© black hole with a = 0.95 and a 10 M© hole with M = 0.01 M© s"^ and 
a — 0.01. These solutions are listed in Table 2, where it can be seen that they have much 
lower neutrino luminosities than 3 Mq black holes with the same values of M, a, and a. 
This is due to the larger size of the more massive holes; this makes the density of the 
accreting gas lower, and thus reduces the pair capture luminosity. 
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5. 



APPROXIMATIONS TO THE ACCRETION RATE 



For those models where the accretion rate is governed by the viscous time scale of the 
disk, it is possible to make a simple estimates of the accretion rate and duration of the event 
provided one knows the radius where the disk forms. This includes most of the models in 
Table 1. One must distinguish however, three classes of solutions: 1) models where the disk 
is being fed at a rate governed by viscous processes outside the region where neutrino losses 
are important - white dwarfs plus black holes and black holes plus helium stellar cores are 
examples; 2) models where the characteristic time scale is the viscous time scale for the disk 
in a region where neutrino losses approximately balance dissipation - merging neutron stars, 
black holes plus neutron stars, and to some extent, the collapsar model; and 3) models 
where the accretion rate is not governed by disk viscosity, but by other circumstances. The 
prime example of this last case is a collapsar where the viscous time in the disk, at least for 
late times and reasonable choices of a, is short compared to the free fall time of matter at 
a much larger radius which falls inward and feeds the disk. The general collapsar model is 
actually a comphcated case that can lie on the boundary of all three classes. 

Consider first the slowest accreting models, those of class 2. The viscous time scale is 
approximately t^sc = t^/i/, where u — aH^flx- Thus we have t^sc = {H / r)~'^Vt~^ . For 
an advection-dominated disk, we have if ~ r, so tyisc ~ W^Q.']^. If we assume that a mass 
Mdisk deposited at a given radius raisk is accreted on the viscous time scale, we can estimate 
the accretion rate as M ~ Mdisk/tvisc- In general, we find for advective disks where neutrino 
losses are negligible 



Assuming the radii for disk formation - approximate values are given by just the Roche 
radius - one obtains the characteristics of merging white dwarf black hole pairs and black 
hole helium core pairs given in Table 1. These agree well with the detailed models (Fryer & 
Woosley 1998; Fryer et al. 1998). The accretion rate spans a large range because high mass 
white dwarfs (and compact helium cores) form their disks at smaller radii than those with 
low mass. 

Next we consider neutrino dominated disks. We assume that neutrino cooling is so 
efficient that any energy dissipated in the disk is quickly radiated away. In other words, 
one assumes that neutrino cooling will produce a thin disk (Shakura & Sunyaev 1973) 
where advection will be unimportant. One can do this by simply substituting the neutrino 
coohng rate (§3) into the standard thin disk equations. There is a question of which of the 



M ~ 0.37aMdiskMl^\2k,9 M© s"^ 



(5-1) 

(5-2) 
(5-3) 
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two types of neutrino cooling described above will be dominant. By solving the thin disk 
equations for each type of cooling in turn, one can then check to see whether this should 
be the dominant cooling term in the resulting disk. Doing so, one finds that the thin disk 
should be quite dense, and, as a result, cooling by pair capture on nucleons should dominate 
over V — V annihilation. The thin disk solution is obtained by equating the cooling rate by 
pair capture to the energy dissipation rate per unit volume {'i/S'n:)MVl\/ H . One then finds: 



T = 1.3 X 10"a°-2Mf°-2i?"°-3 K 


(5-4) 


= 1.7 X 10^«°-^M°-^i?i-3^ cm 


(5-5) 


1.2 X lO^^a-^-^Mf i-^i^-^-^^Mi g cm~3 


(5-6) 


y = 5.6 X lO^a^-^Mf cm s"^ 


(5-7) 


Q = = 2.0 X lO^Mf ^i?-^-^ s"^ 


(5-8) 


Z = 4.4 X 10^^Mii?°-^ cm^ s"^ 


(5-9) 


r = R/V = 2.6 X lO-^a-^-^M^^i?"-^ s 


(5-10) 



where Mi, R, and Mi are the mass of the accreting black hole in solar masses, the radius 
in gravitational radii GMi/c^ (4.43 km for a 3 black hole), and the mass accretion 
rate in solar masses per second. Note that the pair capture rate is proportional to X„uc, 
which is taken to be unity. These equations then provide the estimates given in Table 1 
for the merging neutron stars and neutron star black hole pairs. Again, for an assumed 
disk viscosity a = 0.1 and black hole mass. Mi = 3 M^, they are not far from the detailed 
models (Ruffert & Janka 1998; Eberl, Ruffert, & Janka 1998) provided a disk radius of 50 
km is adopted. This solution is shown in Fig. 1 for comparison with our calculations. 

Finally we consider the collapsar. The angular momentum distribution inside a 
collapsing massive stellar core is unknown. Such calculations as have been done (Heger, 
Woosley, & Langer 1998) indicate an increasing value for / as one moves out in mass, rising 
from a few x 10^^ cm^ s~^ in the inner few solar masses to about 10^'' cm^ s~^ in the outer 
regions. One expects considerable variations from star to star. This range of I corresponds 
to disk radii of 50 - 250 km and viscous time scales (neutrino dominated) of about 0.1 s. 
However, this is considerably shorter than the free fall time for the mass that feeds the disk. 
For a stellar radius that contains most of the mass of, e.g., a 10 helium core, ~ 5 x 10^ 
cm, the free fall time scale 1338 s/p^/^ is several seconds, and this time scale grows longer 
as the accretion continues. The collapsar is probably unable to produce a jet during its first 
5 - 10 s because the infalling matter along the rotational axis sweeps any energy deposited 
into the hole. So a time scale of ~ 10 s seems reasonable for the collapsar. This is consistent 
with detailed 2D models (MacFadyen & Woosley 1998). 
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6. ENERGY CONVERSION EFFICIENCIES 

According to the current paradigm (Rees & Meszaros 1992, 1994, and subsequent 
publications; Katz 1994; Sari & Piran 1997a,b), a gamma-ray burst is formed when 
a relativistic outflow, most likely a jet, suffers collisions both internally and with the 
surrounding medium, producing shocks. Typical values needed for the jet are relativistic 
r's of 100. The total mass and energy of the jet depends upon its opening angle and the 
efficiency for converting its energy to gamma-rays, but should be of order 10~^ — 10~^ 
and lO^'^ — 10^^ erg. Modern gamma-ray burst models thus divide into two categories: those 
that treat the shock interaction and production of gamma-rays and "afterglows" in other 
wavelengths (e.g., Waxman 1997) and those that discuss the energy source itself. We are 
interested here only in the latter problem. 

The main problem in all current models for the GRB energy source is how to 
convert some fraction (albeit small) of the net accretion energy into directed relativistic 
outflow. Two general mechanisms have been proposed - neutrino annihilation and 
magnctohydrodynamical (MHD) mechanisms of various kinds. The former is easier to 
understand and to model and will be dealt with first. 



6.1. Efficiency for neutrino annihilation 

Table 2 shows that a highly variable fraction of the accretion energy will be emitted 
as neutrinos. In general a higher accretion rate and a low viscosity favor high efficiencies, 
LylM(? . A large rotation rate for the black hole also enhances the efficiency. Most of 
these neutrinos come from pair capture on nucleons and are thus solely of the electron 
flavor. The smaller fraction of neutrinos from pair annihilation will be a mixture of three 
flavors. Everywhere in the vicinity of the black hole, but especially along the rotational 
axis, neutrino will encounter neutrino with a large flux. Thus neutrino annihilation, 
i/ -|- P — > e"*" -|- e~, will deposit some fraction of the accretion energy in regions where the 
mass density may be small. These same neutrinos will also drive a wind from the disk by 
their interaction with electrons and nucleons. Since the gravitational potential, neutrino 
temperatures, and fluxes are similar to those for neutron stars, the semi-analytic formulae 
of Qian & Woosley (1996) should hold approximately. This suggests that the wind from the 
disk will be negligible except in the case of neutron star-black hole mergers. There it might 
amount to ~ 0.01 Mq s~^. 

We neglect here the interaction of neutrinos with all other matter and estimate 
the neutrino annihilation efficiency using a simple approximation to the disk geometry. 
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The neutrino luminosity is assumed to be concentrated in the equatorial plane (a good 

approximation given the high temperature sensitivity of the rates) and thus to have a 
luminosity AT:rH{r)q{r) erg cm~^ s~^ where H(r) is the disk scale height and q{r) the 
neutrino luminosity per unit volume evaluated for the temperature and density at radius r. 
The neutrino annihilation at any point above the disk is calculated following the method 
described in Ruffert et al. (1997): 

k "fe fc' "fc' 

(6-1) 

where = cto(Ci + C2)^^^-Jil2n'^c^ml) ^ 1.7 x 10"^^ cmerg-^ s"^ and 
^2 = o'o(C'3)i/ii7i/(47r^c) ~ 1.6 X lO^^^cms"^ are the neutrino cross-section constants 
for electron neutrinos (including geometrical factors). We model the disk as a grid of cells 
in the plane with neutrino mean energy ((e) i^J and luminosity (AL^.). For each pair of 
cells, and every point above the scale height of the disk H{r), we calculate distance {d^) 
from each cell to that point, and the angle (6*) at which the neutrinos and anti-neutrinos 
from the pair of cells interact. The summation over all pairs of cells gives the energy density 
from pair production at that point. 

Each of the neutrino energies and luminosities includes the effects of general relativity 
as the neutrinos travel through the black hole's potential well. Those neutrinos which 
annihilate near the surface of the black hole were emitted further out in the potential well 
and actually gain energy before annihilation. However, the pairs they produce must then 
climb further out of the potential well, an effect that almost exactly cancels out the energy 
gain. In addition, most of the annihilation energy is produced well above the black hole and 
general relativity plays a very minor role (<5%) in the total fireball energy. We have not 
modeled neutrino propagation through a Kerr metric, nor have we modeled the bending of 
geodesies, but it is unlikely that these effects will significantly alter our results. 

For our models, most of the neutrino emission is due to electron capture and the 
luminosity density from pair production (Itpi^i^i)) is dominated by the annihilation of 
electron neutrinos and anti-neutrinos. Although we include the contribution of the n and r 
neutrinos for completeness, they make up less than 10% of the total luminosity density. 

Integrating over the distance above the plane and taking advantage of the cylindrical 
symmetry of our disk: 

27rr / ltA^m)dz (6-2) 
Jo 
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demonstrates the strong focusing of the pair fireball (Fig. 9) For most of the models, over 
half of the energy is injected at equatorial radii ^ 2 x 10^ cm. Further integrating the 
luminosity over the equatorial radius gives the total neutrino annihilation luminosities. 
The net momentum of the pairs produced by this annihilation is directed outward along 
the angular momentum axis, and most of this energy will escape. Table 4 summarizes 
the neutrino annihilation energies and efficiencies for all models studied. The efficiency of 
energy conversion from neutrinos to a electron/positron pair fireball is both a function of 
the total neutrino emission and the distribution of that emission. For example, although the 
M — 0.1 Mq s~^, a — 0.1, a — model has a lower luminosity than its a — 0.01 counterpart, 
because its emission is more centrally concentrated, it converts a larger fraction of the 
neutrinos into electron/positron pairs and its resultant fireball is more energetic. 

Note that the energy conversion efficiencies are extremely high for the M = 10 Mq 
models. For these models, our assumption that the disk is optically thin to neutrinos no 
longer holds (the neutrino mean free path is ~ 1000 — 5000 cm whereas the disk scale height 
is 10^ cm). For the non-rotating model, our estimated luminosity exceeds that of Eberl, 
Ruffert, & Janka (1998) by nearly a factor of 4. Most of this factor arises from the fact that 
our models predict that nearly all of the neutrino energy is produced by electron capture 
(increasing our efficiency by a factor of 3), and hence, is dominated by electron neutrinos 
and anti-neutrinos. Although we have limited our neutrino luminosity by the blackbody 
limit, overestimates in the neutrino temperature no doubt account for the remaining 30-40% 
differences. Our estimates for the M = 10 Mq s~^ models should be taken as an upper 
limit for the neutrino annihilation luminosity, and the actual luminosity could be as much 
as a factor of 5 lower. 

6.2. Efficiencies for MHD energy extraction 

It is also possible to extract energy from either the disk or the black hole by MHD 
processes (e.g., Narayan et al. 1992; Meszaros & Rees 1997; Katz 1997). These are based 
on the expectation that the differential rotation of the disk will rapidly amplify preexisting 
magnetic fields until they approach equipartition with the gas kinetic energy. Examples 
of MHD energy extraction mechanisms include a relativistic wind or jet driven from the 
disk surface (Katz 1997; Meszaros & Rees 1997), Parker instability within the disk leading 
to rcconncction and flares (Narayan et al. 1992), or the Blandford-Znajek mechanism for 
extracting the spin energy of the black hole (Blandford & Znajek 1977). 

Of these, perhaps the one for which it is easiest to estimate the efficiency is the 
Blandford-Znajek (BZ) mechanism. Other MHD processes might have very different 
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efficiencies, particularly since unlike the others the BZ mechanism depends directly on the 
black hole spin, but we give here the BZ-efficiency as a rough estimate. 

All the evolutionary scenarios in Table 1 give black holes which are rotating very 
fast. This is inherent in the fact that they accrete a fraction of their mass from a disk. 
Typical values of the spin parameter a are 0.5 and we adopt that here as a representative 
value. The magnetic field in the accreting plasma is harder to estimate. Here we follow a 
common assumption that the field will rise to some fraction, which we guess might be 10% 
of its equipartition value (i.e., an energy density 1% of pv^). Table 3 gives values of pv^; 
typical values are 10"^° erg cm~^ implying a field strength of 10^"^ — 10^^ gauss. The BZ jet 
luminosity is then 



Entries for our various models are given in Table 6. The estimated luminosities are 
comparable to those for neutrino annihilation listed in Table 5. However, if the magnetic 
field reaches full equipartition with the gas, rather than 1% of pi; ^ as we have assumed, the 
BZ luminosities could be two orders of magnitudes larger than those listed. 

Note that to generate the equipartition field requires not only a large kinetic energy 
density in the vicinity of the of the black hole but that the disk make many revolutions. 
Table 3 gives the number of windings a blob of accreting matter will experience in each 
of our models, which is generally around a~^. The number required to generate the field 
depends upon the primordial field in the accreting matter and the efficiency of instabilities 
in creating radial variations in a magnetic field whose poloidal component is being wound 
by rotation (Katz 1997). We suspect that 10 windings is a gross lower bound. 

7. CONCLUSIONS AND IMPLICATIONS FOR GAMMA-RAY BURSTS 

We have explored accretion at very high rates, 0.01 to 10 Mq into stellar mass black 
holes. The black holes were of both the stationary (Schwarzschild) and rapidly rotating 
(Kerr) varieties. For a disk viscosity a = 0.1 and accretion rates larger than ~ 0.05 M0 s^^, 
we find that a situation of global balanced power comes to exist in the disk interior to 
about 10^ cm where neutrino emission, chiefiy by pair capture on nucleons, approximately 
balances the energy dissipated by viscosity (Table 2). For lower disk viscosity, the necessary 
accretion rate to achieve balanced power is reduced (see a — 0.01, M — 0.01 Mq s"^. Table 
2). Loss of appreciable energy to neutrinos leads to a cooler, denser, and thinner disk than 
the purely advective solution (Fig. 1). The sudden onset of neutrino losses, which are 
very sensitive to the temperature and density, appreciably thins the disk and may help to 
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provide favorable geometry for jet acceleration. 

Temperatures in the inner disk range from 2 MeV to 13 MeV in the models studied 
and densities from 10^ to 10^^ g cm~^ (Table 3). Only for the most extreme accretion 
rates, 1 M© s~^, does the disk start to become optically thick or the neutrino emission 
blackbody limited. Thus one can calculate the neutrino losses and the efficiency for neutrino 
annihilation using a very simple approach. 

We found that the efficiency for converting accreted mass energy, Mc^, into neutrinos 
is highly variable. Greater efficiency is favored by low viscosity and high accretion rate, but 
in no case is the full theoretical limit reached. In some cases the efficiency is very low (e.g., 
0.08% for a = 0.1 and M = 0.01 Mq s'^. The efficiency is also smaller if the black hole 
mass is larger (Table 2). 

We also calculated the efficiency for neutrino annihilation for all of our models and 
found a very large range (Table 4). This range is a consequence of the large fraction 
of accretion energy advected into the hole for accretion rates under 0.05 s~^, which 
reduces the neutrino luminosity, and a reduction in mean neutrino energy for the cooler, 
slower accreting disks, which reduces the cross section for neutrino annihilation. For a 
Schwarzschild black hole, changing the accretion rate by one order of magnitude from 
0.01 Mg s~^ to 0.1 changes the energy deposition by neutrino annihilation by five 

orders of magnitude. Above 0.1 s~^ the efficiency for neutrino annihilation continues 
to increase reaching at least several per cent for accretion rates of 5 s~^ on a rapidly 
rotating black hole. It is impossible to produce a bright GRB (of say 10^^ ergs which, with 
a beaming factor of 100, might resemble a 10^^ erg burst) by neutrino annihilation in any 
model in which the accretion rate is less than a few x 10^^ s^^. On the other hand 
accretion rates over 0.1 s~^, however they may be realized, can give very energetic 
events (Tables 4 and 5). 

With the energy densities calculated for our inner disks we also estimate the efficiency 
for jet acceleration by the Blandford-Znajek (BZ) process using a representative value for 
the spin parameter (0.5) and a guess regarding the the magnetic energy density (1% of pv"^ 
in the inner disk). Typical total GRB energies were around 10^*^ erg (times the beaming 
factor). These estimates are very crude; in particular, the field energy density in the disk 
could be 1%, 100%, or 0% of equipartition. We give these estimates only to show that the 
kinetic energy density in the disk is sufficient to anchor an adequate field to power a bright 
(beamed) GRB. 

However, we note that these large fields might require many orbits of the disk to 
generate. The number of disk windings is small unless the viscosity is initially low (Table 
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3). One can envision a disk that has essentially zero viscosity until the field is wound up 
to a significant fraction of equipartition. Then the estimates of field strength in the inner 
disk (Table 6) would be higher. We see two solutions to the GRB problem emerging - 
high viscosity disks that make bursts by neutrino annihilation, and low viscosity disks that 
extract rotational energy from the black hole by MHD processes. Our understanding of the 
disk physics is inadequate to distinguish these at the present time. We have concentrated 
on neutrino annihilation here simply because it is easier to calculate. 

The implications of our work for GRB models arc best discussed on a case by case 
basis (sec Tables 5 and 6). The event rates in Table 5 arc representative and are taken from 
the forthcoming paper by Fryer, Woosley, & Hartmann (1998). 

7.1. Merging Neutron Stars and Black Hole Neutron Star Pairs 

Our calculations here qualitatively confirm the more accurate, but specialized work of 
Ruffert & Janka (1998) and Eberl, Ruffert, & Janka (1998). Given that a disk forms at ~50 
km, eq. 5.9 gives the correct approximate accretion time scale, 0.1 s, provided a = 0.1 and 
the black hole mass is 3 . For other values of a and black hole mass, our results suggest 
how their results might be scaled. Our neutrino luminosity, about 10^^ erg s~^, neutrino 
annihilation rate, 5 x 10^° erg s~^, and total energy available for the burst, 5 x 10^^ erg, 
agree with the detailed calculations of merging neutron stars. 

For neutron star plus black hole mergers, our calculations suffer from the disk becoming 
optically thick. The time scale for accretion, still 0.1 s, agrees well, and for the larger disk 
mass (Table 1), it is clear that greater energies and efficiencies will be developed. The 
entries in Table 5 come from logarithmic interpolation in Table 4 to obtain the neutrino 
annihilation luminosity for an accretion rate of 5 s~^, then arbitrary division by 4 to 
correct for the decrease in neutrino luminosity and cooling of neutrino energy one expects 
for a (optically) thick disk. These numbers are then consistent with Eberl, Ruffert, & Janka 
(1998). 

The chief new result of our calculations here is an estimate of the effect of using Kerr 
geometry. The 3D numerical calculations of Eberl, Ruffert & Janka (1998) were all carried 
out for a non-rotating black hole. Here we see that a Kerr parameter of only 0.5 increases 
the efficiency of neutrino (emission and) annihilation by a factor of 4 to 6. This degree of 
rotation is reasonable for a black hole made by merging neutron stars, and, depending on 
its formation process, might be reasonable for a black hole merging with a neutron star. 

Of all the models in Table 5, only the merging neutron stars and black hole-neutron 
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star pairs have a characteristic time scale much less than one second. They are thus the 
only model capable of explaining the abundant subclass of short (average duration 0.3 s) 
hard GRBs (Fishman & Meegan 1995). We have no information on the optical counterparts 
of these short bursts, and do not know if they are associated with star formation. 

Because the initial conditions are so uniform, especially for merging neutron stars, one 
might expect that the jet has a nearly imique energy and duration. GRB diversity in this 
case would solely reflect the variable characteristics of the interstellar medium in which the 
burst was embedded. 



7.2. Collapsars 

Here our results are more sensitive to uncertain parameters and a diversity of outcomes 
is possible. In addition to the uncertain disk viscosity which plagues all accretion powered 
models, the accretion rate is particularly sensitive to the angular momentum distribution 
of the accreting star. This is poorly known and likely to be highly variable both within a 
star and from star to star. The radius where the disk forms varies as (eq. 5.8) and the 
accretion rate as an even higher power. Fortunately, if the disk viscosity is adequately high, 
and if a steady state disk forms, the accretion rate may be determined by simpler physics 
- the rate at which matter is delivered to the accretion disk by the collapsing star. Our 
estimate comes from taking the free fall time of the star divided into that fraction of the 
star's mass which forms a bound disk. We estimate 10 s for the former and perhaps a few 
solar masses for the latter. The frequency of events is also very uncertain. In Table 5 we 
have taken from 0.5% to 5% of the supernova rate. 

The formation of a bound disk requires dissipation. The disk of a coUapsar may be 
unique in nature in that it is assembled from optically thick matter falling in, essentially 
freely from infinity. In the absence of neutrino losses and photodisintegration (which only 
occur inside a few hundred km), the net binding energy of the disk must be carried away 
by outfiowing matter. This can be a considerable amount of energy. For example, a 1 
Mgdisk forming at 250km (/ ~ 10^^ cm^ s~^) must release a binding energy of 3 x 10^^ erg. 
Half of this is rotational energy and most of the rest is internal energy. But if even a 
small fraction goes into outflows, they could power an energetic supernova like SN 1998bw 
(Woosley, Eastman, & Schmidt 1998). Because a lot of the matter that flows out will have 
been at small radii and very hot, its composition may be rich in ^^Ni. Our simple steady 
state calculations here cannot model the multi-dimensional physics of disk formation and 
bi-directional flow. Two dimensional studies are in progress (MacFadyen & Woosley 1998). 
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But provided enough dissipation occurs, by mass ejection, viscous interaction, 
photodisintegration, and neutrino losses, something resembhng a steady state accretion 
disk will form. As the accretion proceeds the black hole will be spun up. Our results show 
(Table 4) that the energy available for jet production from neutrino annihilation is highly 
variable ranging from essentially zero, if the accretion rate is 0.01 M^^^ s^^ or less, to over 
Iq52 j£ black hole rotates rapidly and the accretion rate exceeds 0.1 s~^ for 10 
s (Table 5). This makes coUapsars potentially the most energetic and most frequent of all 
GRB models. If the beaming factor is ~100, as preliminary calculations by MacFadyen & 
Woosley (1998) indicate, GRBs of equivalent energy up to 10^^ erg might be explained. 

It should be noted that in this and all models with long time scales, the energy available 
for a GRB is the GRB duration times the jet luminosity - not necessarily the total energies 
in Tables 5 and 6. This is another reason why only merging neutron stars and black holes 
can make energetic GRBs shorter than 1 s. 

7.3. White Dwarfs Plus Black Holes 

For these models the accretion rate is expected to be lower. The Roche radius gives a 
range of disk sizes and, for a given disk viscosity and white dwarf mass, this sets the range 
of accretion rate (Table 1). For the lower values, the energy from neutrino annihilation 
is negligible. For the higher mass white dwarfs which give accretion rates of 0.07 s~^ 
for 15 s, the energy is dramatically greater. If the black hole rotates, the yield is further 
increased to perhaps 10^*^ erg. We expect the beaming factor to be less here than in the 
collapsar, but bursts of up to 10^^ erg might be explained, more if the Blandford-Znajek 
mechanism is effective (Table 6). However, unless these events occur more frequently than 
black hole-neutron star mergers, they are not likely to be the leading cause of GRBs. 

7.4. Black Holes and Common Envelope Evolution 

As with the collapsar which it resembles, this is a model where critical parameters 
remain poorly determined. The three dimensional evolution of a black hole merging with 
a helium core (or helium star) of comparable or greater mass has not been studied. One 
expects a great deal of the accretion to occur along the poles. The radius at which the 
hehum disk forms is poorly known. If it is as small as 10^ km for at least 1 of helium, 
accretion rates of 0.1 s~^ may be reahzed. If the black hole further has a Kerr parameter 
of 0.5, the neutrino annihilation energy would be ~ 10^^ erg. If the rotation of the black 
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hole becomes even faster as a consequence of the merger, the energy could be increased still 
more. Because of the large mass and angular momentum in the merger, the channeling of 
the jet may result in tighter beaming, comparable to the coUapsar and more focused than 
merging compact objects. 
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TABLE 1 



Model 


Maccrete 


Radius 


M 


j 


Duration 


Gravity 






(km) 


(Mjs) 


(10^6 cm2 s-i) 


(sec) 


Waves 


n*+n* 


0.1 


50 


1 


4 


0.1 


Yes 


n*+BH 


0.5 


50 


5 


4 


0.1 


Yes 


CoUapsar 


2 


50 - 250 


0.1 


5 - 10 


10-20 


No 


BH+WD 


1 


(1-5) X 10^ 


0.01 - 0.07 


50 - 150 


15 - 150 


No 


BH+He core 


2 


(1- 10)xlO^ 


0.01 - 0.1 


50 - 200 


15 - 500 


No 



Masses are for accretion through a disk. The total accretion rate, e.g. in the coUapsar and 
He core models, is greater because of mass infall along the poles. The assumed mass of the 
black hole is 3 M in all cases and the disk viscosity, a = 0.1. 
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TABLE 2 





Model 








Luminosities (10^^ 


erg/ s) 






M 


a 


a 


M 








Lv,URCA 


Lphoto 


^diss 


Lent 


0.01 


0.1 





3 


0.015 


8.6e-4 


3.4e-3 


0.012 


0.200 


0.859 


0.644 


0.01 


0.03 





3 


0.089 


5.0e-3 


0.012 


0.078 


0.200 


0.738 


0.450 


0.01 


0.01 





3 


0.650 


0.036 


0.012 


0.638 


0.200 


0.724 


-0.124 


0.01 


0.1 


0.5 


3 


0.036 


2.0e-3 


9.9e-3 


0.026 


0.200 


1.24 


1.01 


0.01 


0.01 





10 


0.049 


2.7e-3 


9.3e-3 


0.040 


0.200 


0.431 


0.181 


0.05 


0.1 


0.5 


3 


1.65 


0.018 


0.36 


1.34 


1.00 


6.48 


3.86 


0.1 


0.1 





3 


3.35 


0.019 


0.39 


2.96 


2.00 


9.18 


3.88 


0.1 


0.03 





3 


6.96 


0.039 


0.076 


6.89 


1.99 


7.94 


-1.01 


0.1 


0.01 





3 


6.15 


0.034 


7.7e-3 


6.14 


2.00 


7.44 


-0.71 


0.1 


0.1 


0.5 


3 


8.03 


0.045 


1.03 


7.00 


2.00 


12.96 


3.02 


0.1 


0.1 


0.95 


3 


46.4 


0.26 


2.61 


43.8 


1.99 


36.8 


-11.5 


0.1 


0.1 


0.95 


6 


26.2 


0.15 


7.24 


19.0 


2.00 


31.1 


3.20 


1.0 


0.1 





3 


86.3 


0.048 


0.66 


85.7 


19.9 


90.9 


-14.9 


1.0 


0.1 


0.5 


3 


142. 


0.078 


1.53 


140. 


19.9 


142. 


-19.6 


10.0 


0.1 





3 


781. 


0.043 


53.3 


728. 


200. 


836. 


-135. 


10.0 


0.1 


0.5 


3 


1285. 


0.071 


174. 


1111. 


200. 


1263. 


-211. 
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TABLE 3 





Model 






Values at horizon 




M 


a 


a 


M 


P 


1 


pV 




iVi0 b 






iVi0 


1 0-^^ rr rim 
lU g, Clli 


lO^o K 






0.01 


0.1 





3 


0.060 


3.90 


0.53 


7.1 


0.01 


0.03 





3 


0.084 


3.66 


0.75 


23.4 


0.01 


0.01 





3 


0.23 


1.91 


2.04 


90.2 


0.01 


0.1 


0.5 


3 


0.080 


4.61 


0.70 


7.6 


0.01 


0.01 





10 


9.32e-3 


1.96 


0.083 


19.7 


0.05 


0.1 


0.5 


3 


0.43 


6.67 


3.84 


9.8 


0.1 


0.1 





3 


0.69 


6.46 


6.14 


9.6 


0.1 


0.03 





3 


1.67 


3.47 


15.0 


31.2 


0.1 


0.01 





3 


2.10 


2.04 


18.8 


87.9 


0.1 


0.1 


0.5 


3 


1.07 


7.42 


9.45 


10.2 


0.1 


0.1 


0.95 


3 


12.0 


5.83 


98.7 


14.0 


0.1 


0.1 


0.95 


6 


0.88 


9.04 


7.46 


10.2 


1.0 


0.1 





3 


12.6 


5.71 


113. 


9.9 


1.0 


0.1 


0.5 


3 


20.1 


6.35 


178. 


11.5 


10.0 


0.1 





3 


86.6 


11.4 


774. 


6.6 


10.0 


0.1 


0.5 


3 


135. 


13.1 


1197. 


7.4 
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TABLE 1 



M 


a 


a 


M 


U 




efficiency 








Mo 


(10^1 erg s-^) 


(10^1 erg s-^) 


(%) 


0.01 


0.1 





3 


0.015 


3.9 X 10 


0.0003 


0.01 


0.03 





3 


0.089 






0.01 


0.01 





3 


0.650 


9.0 X 10-6 


0.001 


0.01 


0.1 


0.5 


3 


0.036 


5.9 X 10"^ 


0.002 


0.01 


0.01 





10 


0.049 


6.4 X lO"'^ 


10-^ 


0.05 


0.1 


0.5 


3 


1.65 


1.8 X 10-3 


0.11 


0.1 


0.1 





3 


3.35 


3.0 X 10-3 


0.09 


0.1 


0.03 





3 


6.96 






0.1 


0.01 





3 


6.15 


8.0 X 10-^ 


0.01 


0.1 


0.1 


0.5 


3 


8.03 


0.039 


0.5 


0.1 


0.1 


0.95 


3 


46.4 


2.0 


4.2 


0.1 


0.1 


0.95 


6 


26.2 


0.79 


3.0 


1.0 


0.1 





3 


86.3 


0.56 


0.6 


1.0 


0.1 


0.5 


3 


142 


3.5 


2.5 


10.0" 


0.1 





3 


(781) 


(200) 


(26) 


10.0" 


0.1 


0.5 


3 


(1280) 


(820) 


(64) 



The assumption tliat tlie neutrinos are optically thin breaks down for accretion rates of 10 
MgS-^ and above. The neutrino annihilation luminosities and energies listed for these high 
accretion simulations are upper limits. 



-34- 



TABLE 5 



Model 


Duration 


a 


log Kv 


log E^p 


Rate 




(s) 




(erg s-^) 


(erg) 


(Myr-i)" 


n*+n* 


0.1 





50.7 


49.7 


0.1-5 






0.5 


51.5 


50.5 




n*+BH 


0.1 





(52)^ 


(51)" 


0.1-50 






0.5 


(52.6)^ 


(51.6)'' 




CoUapsar 


10 





48.5 


49.5 


^ 2000 






0.95 


51.3 


52.3 




BH+WD 


15-150 





43.6 - 47.7 


46-49 


0.1-50 






0.5 


44.8 - 48.8 


47-50 




BH+He core 


15-500 





43.6 - 48.5 


46-50 


1-1000 






0.5 


44.8 - 49.6 


47-51 





"This rate assumes a supernova rate of 0.02 yr 
^Rough estimate owing to optically thick disk. 
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TABLE 6 

Model" Duration B15 log Lbz log Ebz 

(s) (10"^^ gauss) (erg s^-*^) (erg) 

n*+n* OA 10 51 50 

n*+BH 0.1 10 51 50 

CoUapsar 10 1 49 51 

BH+WD 15-150 1 49 50 

BH+Hecore 15-500 0.3 48 51 



"Assuming a=0.5 in all cases 



-36- 



Fig. 1. — A neutrino-cooled accretion disk solution (solid line) for an accretion rate 
M = 1 Mq s~^, black hole spin a = 0, viscosity parameter a — 0.1, and black hole mass 
M = 3 Mq. The six panels show the density p, temperature T, radial velocity V, disk 
height H, angular velocity Q, and specific angular momentum /. Also shown for comparison 
purposes are a pure advection solution with no cooling (dotted line) and an analytic thin disk 
solution where URCA-process neutrino cooling balances viscous dissipation (dashed line). 

Fig. 2. — Additional parameters for the solution shown in Fig. 1. Note that here the dashed 
and dotted lines have different meanings than in Fig. 1. The upper left panel shows the 
cooling rate due to neutrinos and photodisintcgration (solid line), the heating rate due to 
viscous dissipation (dotted line), and the rate at which energy is added to the entropy of 
the gas (dashed line). The upper right panel shows / = 1 — Qcooi/Qdiss', the feature around 
log r — 6 occurs because Qdiss briefly becomes negative here. The middle left panel shows the 
neutrino luminosity per unit radius. The middle right panel shows the cooling and accretion 
time scales. The two bottom panels show the pressure components: total pressure (solid), 
gas pressure (dotted), pressure due to radiation and pairs (short-dashed), and degeneracy 
pressure (long-dashed). 

Fig. 3. — Solutions for four accretion rates: M — 0.01 (short-dashed), 0.1 (solid), 1 (dotted), 
and 10 Mq s"-*^ (long-dashed). All four have a — 0, a — 0.1, and M — 3 Mq. The panels are 
as in Fig. 1, except for the two bottom panels which show the neutrino luminosity per unit 
radius and /. 

Fig. 4. — Solutions for three values of the viscosity parameter a = 0.1 (solid), 0.03 (dotted) 
and 0.01 (dashed). All three solutions have M = 0.1 Mq s~^, a = 0, and m — 3 Mq. 

Fig. 5. — Same as Fig. 4, except that here the three solutions have M — 0.01 Mq s~^. 

Fig. 6. — Solutions for three values of the black hole spin: a — (solid), 0.5 (dashed), and 
0.95 (dotted). All three solutions have M = 0.1 Mq s^^ a = 0.1, and M = 3 Mq. 

Fig. 7. — Same as Fig. 3, except that here the four solutions all have a = 0.5. 

Fig. 8. — The Kerr parameter a as a function of the fractional increase in the mass of a 
black hole accreting from a thin disk. 

Fig. 9. — (a) Pair annihilation luminosity per cm versus radius for M = 0.1 Mq s~^ and 
Mbh = 3.0 Mq. In addition, the total integrated luminosity out to each radius is plotted 
for the same models. These models are representative of the entire set, for all of which over 
half of the energy is injected in a narrow beam with equatorial radii ^ 2 x 10^ cm. (b) Pair 
annihilation luminosity per cm versus distance along the angular momentum axis. Again, 
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the total integrated luminosity out to each scale height above the plane is plotted for the 
same models. 



-39- 




6789 6789 

log r log r 



-40- 




-41 - 




-42 - 




6 7 8 

log r (cm) 



9 



6 7 8 

log r (cm) 



9 



-43- 




6 7 8 

log r (cm) 



9 



6 7 8 

log r (cm) 



9 



-44- 




-46- 




Total Luminosity (erg s~') 



I L I — I — I — I 1 r- 



O O 
II II 
CO (0 



^ ^ O 

odd 
II II II 
« is a 



_I lli^ I I I I I I L 




o o o 

(,_s ,_uio Sja) zp (<i'z'j)-ii/jug ruojiBnuiuuv JIBJ 



